Evolution of a black hole at the center of GRB 
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Abstract 

Using a simplified model of a black hole-accretion disk system which is dominated by Poynting 
flux, the evolution of the central black hole which is supposed to be powering GRB is discussed. 
It is demonstrated explicitly that there is a lower limit on the angular momentum parameter for 
a given GRB energy. It is found that the most energetic GRBs can only accommodate relatively 
rapid-rotating black holes at the center. For a set of GRBs for which the isotropic energies and 
Tggs are known, the effect of the disk mass and the magnetic field on the horizon are discussed 
quantitatively. It is found that the magnetic field has little influence on the energy but affects the 
GRB duration as expected. The role of the disk mass is found to be significant in determining 
both the energy and the duration. 
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I. INTRODUCTION 



The discovery and the recent observations followed afterwards indicate that the central 
region powering gamma ray burst is rather compact in size, less than 10^ cm, from the 
studies of energies and temporal structures |1|]. Merging binary compact objects, hypernovae 
and rotating black holes have been considered to be among the viable candidates for a GRB 
central engine. The common feature of these models is the formation or the existence of a 
black hole at the central region. However, the gamma ray bursts or the afterglows followed 
do not provide any direct observational informations on the central region and we do not 
know the exact nature of the black hole at the center. Hence although it is very important 
and interesting to investigate the central object itself, we should rely on the indirect method 
of using a particular model, in which the physical properties of black holes can be inferred 
from the observational data. 

Among the models proposed so far, we choose a model in which the rotational energy 
of the black hole is responsible for powering GRBp|. The mechanism of tapping rotational 
energy from the black hole has been known as Blandford-Znajek mechanism0[^, in which 
the rotational energy is extracted out to the loading region via the magnetic flux which 
threads on the horizon. It is easy to demonstrate that the strong magnetic field of ~ 
lO^^gauss is consistent with the essential features of GRBs|0. 

It is well known that the black hole itself cannot keep the magnetic field on the horizon. 
It disappears very rapidly if not the environment keeps them from disappearing. Recently, 
it is shown that the magnetic flux on the strongly magnetized object can be maintained 
during the collapsing process together with the electric charge onto the black hole[^. When 
the gravitationally unstable object collapses into black hole the most natural environment 
is the accretion disk/torus which emerges together with the central black hole and it can 
provide the magnetoshpere which keeps the magnetic flux from disappearing. Because of 
the accretion the energy and the angular momentum are carried into the black hole while 
the Poynting flux carries away energy and angular momentum out of the black hole[^], the 
evolution of a black hole depends not only on the Poynting flux but also on the accretion. 

In this work we make an attempt to infer the evolution of the central black hole during 
the gamma rays are bursting using a simple model of a black hole-accretion disk system 
suggested in 0, in which the accretion is dominated by the Poynting flux0[^]. In section 
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II, the simple minded model is sketched with emphasis on the parameters which govern the 
evolution of the black hole. We take the initial values of the magnetic field on the horizon 
and the disk mass as parameters while the black hole mass is initially taken to be a typical 
mass, 7M0, suggested from the black hole binary systems. The evolution of the angular 
momentum parameter a is discussed in detail in section III. With finite mass of accretion 
disk, Md < 7Mq, the initially rapid- rotating black hole is expected to remain rotating 
even at the end of GRB. But the initially slow-rotating black hole is found to eventually 
stop rotating with GRB. In section IV, the total energy out of the rotating black hole is 
calculated and compared to GRBs for which the isotropic energies are known |1T0|. The lower 
limits of the angular momentum parameter of the black holes for the corresponding GRBs 
are calculated. For the most energetic GRBs, Eiso > -IMq, the lower limit is found to be 
rather high, a(0) > 0.3. In section V, the evolution of the black hole is discussed for GRBs 
with known Tqq. Identifying Eiso as 90% of energy extraction and T90 as the time taken to 
extract out Eiso, the corresponding sets of Bh{0) and MoiO) are determined to discuss the 
evolution of the corresponding black holes. The discussions are given in section VI. 



II. BRIEF SKETCH OF A MODEL 



In this work, the energy of GRB is supposed to be powered by Poynting flux out of the 
rotating black hole surrounded by the magnetized accretion disk0]. Fig. 1. To see the 
effects of the magnetic field on the evolution of the black-hole-accretion disk, it is assumed 
that the accretion flow is driven mainly by the strong magnetic field which carries angular 
momentum by Poynting fiux0||^. Then we can have an analytic formulation of the black 
hole evolution. The essential feature of the evolution model discussed in ref |^ is briefly 
reviewed in the following. 

The evolution rates of the black hole mass (M) and the angular momentum (J) are 
determined both by the energy and the angular momentum accreted, which increase the 
mass and the angular momentum, and also by the Blandford-Znajek power in the opposite 
direction. From the energy and angular momentum conservation, we get the rates of change 
for the mass and the angular momentum given by 

M = -Pbz + M+E, (1) 
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j = -^+M+l (2) 

where M+ is the mass accretion rate and the specific energy and angular momentum of the 
accreting matter are denoted as E and I respectively. Since Pbz is proportional to M_|_ in 
this particular model[0|, one can get 

Pbz = {E - Pe)M+, j = MM+{Zi - Pi), (3) 



where 



-2Z + a^ ~ _ f{h)H^ 

I =1 Pe — 



E = , = + iaiZfil + 2/Z)], (4) 



, Z^-2~a^fZ+~a'' 
l = MZi{a), Zi{a) = (5) 

/Z(Z2 -3Z + 2~a^) 



Pi=^-^^[l + {h/Z)\l + 2/Z)l (6) 



and 



/(/i) = — -)arctan/i- 1], h = —, H = l + ^l-dfi (7) 
h h H 



It is assumed that the inner edge of the accretion disk is the last stable orbit |]TT| defined by 



Tin = ZM, (8) 

Since the presence of an accretion disk with an appreciable magnetic field is essential for 
the magnetic field on the horizon, the evolution of the accretion disk is also responsible for 
the evolution of the magnetic field on the disk and Bh- 

For the numerical calculations in this work, we assume that the effects of mass loss from 
the disk can be incorporated into the the time dependence of the magnetic field in the 
following form: 

5| = BUO)D{t), (9) 

where 

D{t) = l-{f M+)/Mj,{Q). (10) 
Jo 

We take D{t) to be vanishing as the total accreting mass becomes the initial disk mass, 
which is supposed be typically of solar mass for an accretion disk emerging out of binary - 
merging processes. 



III. EVOLUTION OF ANGULAR MOMENTUM 



The rotation of a black hole can be described either by angular momentum J, or specific 
angular momentum a(= J /M) or angular momentum parameter a(= a/M). In this work we 
will take the angular momentum parameter a which is a dimensionless quantity to represent 
the rotation of the black hole. 

The evolution of the angular momentum parameter a can be obtained from the angular 
momentum evolution as 

j 
1 

Then, using Eq. (|), we get 

a 

where the sign factor A defined by 

A={Zi~ Pi) - 2a{E - Pe) (13) 

determines whether a is decreasing or increasing. As shown in Fig. 2, the sign factor is 
negative for the entire range of a. It means that the rate of change of the angular momentum 
parameter a, is always negative in this model described in the previous section. It should be 
noted that this sign factor is independent of the parameters of the model. It is the intrinsic 
feature of the model. 

The sign factor has a maximum at a ^ = .57 but still less than zero. It is analogous 
to the energy barrier to be overcome in particle dynamics. If the evolution starts with 
a > .57 one should wait long time until the black hole stops rotating. It is because when it 
is crossing the maximum value, (im, the process becomes very slow. With the finite mass of 
the accretion disk, the evolution ends when the disk disappears completely into the black 
hole. In this case however the evolution may end up with the rotating black hole as one can 
see in Fig. 3. On the other hand, if the evolution starts with a < .57 the slow-down process 
accelerates in the beginning and the black hole eventually stops rotating in a much shorter 
time. However it depends also on the evolution of the disk. Since the life time of the disk 
depends essentially on the initial mass, we can define, for a given mass of accretion disk, 
ac(Fig. 4): a black hole start with d < dc ends up with non-rotating black hole. 
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IV. ENERGY OF A GAMMA RAY BURST 



In this work, the energy of GRB is supposed to be powered by the Poynting flux out 
of the rotating black hole at the center of a black hole - accretion disk system. Essentially 
the energy is a part of the rotational energy of the black hole. The rotational energy of the 
disk is accreted into the black hole and then it is processed to be a part of the black hole's 
rotational energy. Therefore the source of the energy delivered to GRB is not only the black 
hole's initial rotational energy but also the rotational energy of the disk. 

The duration of the GRB is identified to be the duration of Poynting flux from the system. 
The Blandford-Znajek process becomes ineffective either when there is no magnetic field on 
the black hole or it does not rotate. If the system starts with the larger disk mass, the hfe 
time is longer for a given initial magnetic field. Therefore the energy carried out until the 
Blandford-Znajek process stops is expected to be larger for the larger disk mass system. 
Numerical calculations in Fig. 5 shows that the energy is increasing with the initial disk 
mass, Md{0). 

However for d < dc, the energy reaches maximum at a certain value of disk mass, M^, 
and shows httle change for the larger disk mass. It is because the rotation of the black hole 
slows sufficiently down for the Blandford-Znajek process to be ineffective as discussed in the 
previous section. When the system starts with Mu{0) > the black hole is slow down 
substantially before the life time of the accretion disk. And the GRB duration is determined 
by the time when black hole stops rotating not by the life time of the disk. The radial current 
in the disk vanishes and there is no magnetic breaking hence no accretion onto the black 
hole. For d > d^ the disk disappears into the black hole while the remaining black hole is 
still rotating. But there is no Poynting flux because the magnetic field supported by the 
disk disappears altogether. Then the duration of GRB is determined by the disk life time. 

Since we are interested in the system of black hole - accretion disk which emerges during 
a collapse or merging process of stellar objects, Md might not be much greater than the 
black hole mass in the system. And for numerical calculation we take the upper limit to be 
7Mq in this work. It gives upper limits of the energy out of the system as shown in Fig. 6 
for given initial values of the angular momentum parameter. 

Compared to the isotropic energies of GRB, one can estimate the lower limits of the 
initial angular momentum parameters, di, of the black holes at the center of the GRB in 
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this model, Table I. One can see that the energetic GRBs require rapidly rotating black 
holes at the center: for example a(0) > 0.45 for GRB 990123. 

V. EVOLUTION OF A BLACK HOLE 

The effect of the magnetic field on the total energy out of the system is not significant as 
shown in Fig. 7. But the strength of the magnetic field is directly involved in determining 
the rate of energy carried out along the magnetic flux. The stronger magnetic field extracts 
energy more rapidly than the weaker magnetic field. Hence the detailed variation of the 
GRB duration is expected to be related to the magnetic field structure of the system. In 
this work it is parameterized by the initial value on the horizon, Bh{0)- 

The determination of the duration of GRB from the observational data is relatively 
complicated than the isotropic energy. In this work, as a first trial, we take Tgo as the 
duration for 90 % of the total Poynting flux energy to be carried out. Among the GRBs 
with inferred isotropic energies in Table I, we choose six of them in Table II for which Tgos 
are well defined. 

The evolution of a black hole is determined by the mass and angular momentum. For 
the system with a(0) = .8, it is found that the final mass of the black holes are found to be 
increasing: AMbh = •23Mq upto 2.8Mq. Since a(0) = .8 > dm, the final black holes are 
expected to be rotating even after GRB with smaller angular momentum parameters than 
the initial values. For a(0) = .3, which is smaller than the lower limits of GRB 990123, 
990506, 991216, 000131 which require more rapidly rotating black holes at the center, the 
less energetic GRBs, 990510 and 991208, can accommodate the slowly rotating black holes. 
In this case, however, the appropriate disk masses are found to be smaller than Mf^ and 
the black holes remain rotating after GRB. It is because a(0) is greater than dc for the 
corresponding disk mass. 

VI. DISCUSSION 

Using a schematic model for a central engine of GRB, which consists of a black hole and 
an accretion disk0] at the center, the evolution of the central black hole is discussed. The 
accretion is assumed to be dominated by the Poynting flux out of the disk and the GRB 
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is supposed to be powered by the Poynting flux which extracts out a part of the rotational 
energy of a black hole. It is found that the evolution of the rotation parameterized by a{t) 
shows different patterns depending on the initial value, a(0). For a given disk mass, the 
black hole with a(0) < .57 much more rapidly approaches to non- rotating black hole than 
with a(0) > .57. It is also shown that there is a maximum energy for GRB for a given 5(0). 
Hence one can infer the lower limit of the angular momentum parameter for the central 
black hole. The most energetic GRBs, are found to be able to accommodate only rapidly 
rotating black holes , a(0) > 0.4. 

The effect of the magnetic fleld on the total energy for GRB is found to be not signiflcant 
compared to the disk mass. However since the stronger magnetic fleld extract energy more 
rapidly than the weaker magnetic fleld, the detailed variation of the GRB duration is found 
to be due to the magnetic fleld structure of the system. The role of the disk mass in this 
model is signiflcant both in determining the energy and the duration of GRB. Within the 
range of the parameters used to flt a set of GRBs for which isotropic energy and Tgo are well 
determined, the flnal black holes are found to become more massive than the initial values but 
with smaller angular momentum parameters. This observation is consistent with the general 
feature expected for a system of a black hole - accretion disk with flnite size and life time. 
As it is mentioned the analysis is limited to GRBs with Tgo determined in this work. But 
provided with a systematic way of determining GRB duration time from the observational 
data, statistically the more meaningful conclusion can be made on the evolution of the black 
hole at the center of GRB and it remains as a future work. The characteristic feature of 
this model is the sign factor A determined by the magnetically dominated accretion disk. 
However it depends on the several simpliflcations which are subject to be verifled. For 
example, we use the relation of the fleld components suggested by Blandford which requires 
a justiflcation if it can be used in the relativistic formulation especially for a rapidly rotating 
black hole. Also the identiflcation of the angular velocity of the fleld line as the Keplerian 
angular velocity in this model also needs a valid justiflcation p. 

This work was supported by Hanyang University, Korea made in the program year of 
2001. 
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GRB 


z 


di [cm] 


Eiso[erg] 


ai 


970228 


0.695 


1.403C+28 


1.42C+52 


0.14 


970508 


0.835 


1.757e+28 


5.46e+51 


0.09 


970828 


0.958 


2.082e+28 


2.20e+53 


0.31 


971214 


3.418 


9.877e+28 


2.11e+53 


0.31 


980613 


1.0969 


2.459e+28 


5.67e+51 


0.09 


980703 


0.9662 


2.103e+28 


6.01e+52 


0.21 


990123 


1.6004 


3.925e+28 


1.44e+54 


0.45 


990506 


1.30 


3.037e+28 


8.54e+53 


0.42 


990510 


1.619 


3.982e+28 


1.76e+53 


0.30 


990705 


0.84 


1.770e+28 


2.70e+53 


0.33 


990712 


0.433 


7.927e+27 


5.27e+51 


0.09 


991208 


0.707 


1.429e+28 


1.47e+53 


0.28 


991216 


1.02 


2.250e+28 


5.35e+53 


0.39 


000131 


4.500 


1.369e+29 


1.16e+54 


0.44 


000301C 


2.034 


5.269e+28 


4.64e+52 


0.19 


000418 


1.118 


2.523e+28 


8.29e+52 


0.24 


000926 


2.037 


5.280e+28 


2.97e+53 


0.34 



TABLE I: The isotropic energies of GRB [10] and the lower limits of angular momentum parameter, 
ai. 
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GRB 


990123 


990506 


990510 


991208 


991216 


000131 


red shift, z 


1.6004 


1.3 


1.619 


0.707 


1.02 


4.5 




0.7167 


0.395 


0.0678 


0.0644 


0.2417 


0.572 




24.34 


57.10 


25.80 


39.84 


7.51 


9.09 


a(0) = 1 


Bu{Qi)\lQi^^aauss\ 


5.19 


2.56 


1.49 


1.17 


5.44 


7.69 


Mn(0)\Mn] 


3.28 


1.54 


0.16 


0.15 


0.82 


2.47 




8.71 


7.71 


7.04 


7.03 


7.33 


8.24 




0.80 


0.87 


0.98 


0.98 


0.92 


0.83 


a(0) = 0.8 


-Bii(0)[10"^\9r/(/.s.s-] 


5. 96 


3.21 


2.21 


1.73 


7.29 


9.13 


Mb(O)[M0] 


4.36 


2.34 


0.39 


0.37 


1.41 


3.44 




9.75 


8.46 


7.24 


7.23 


7.87 


9.16 


af 


0.71 


0.74 


0.79 


0.79 


0.76 


0.72 


a(0) = 0.3 


BHiO)[10^^ gauss] 






7.31 


5.68 






Mb(O)[M0] 






0.59 


0.56 












7.42 


7.40 






af 






0.22 


0.23 







TABLE II: The initial disk mass and the magnetic field of black hole for GRBs for which Tgos are 
well defined. and are the final black hole mass and the angular momentum parameter 

respectively. The energetic GRBs can not accomodate the black hole with a(0) = 0.3. 
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GRB 




Black Hole Disk 



FIG. 1: Skematic representation of the black hole - accretion disk system powering GRB(only 
upper half plane is shown). The rotational energy extracted out by poynting flux Pbz is supposed 
to power GRB. The evolution of the black hole with the mass M and the angular momentum J is 
determined by Pbz and the accretion rate M+. 
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FIG. 3: The evolutions of the angular momentum parameter in time with M(0) = 7Mq and 
Md(0) = 3Mq. 
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FIG. 4: Oc for a given mass of accretion disk in the unit of solar mass M solar- 



15 



2.5 



2 - 



1.5 



0.5 



a(0)=0.1 

a(0)=0.2 

— - a(0)=0.3 

-- a(0)=0.4 

- a(0)=0.5 

- - a(0)=0.6 

- - a(0)=0.7 

- -- a(0)=0.8 



/// y 



/ // y^ 



/ , y y 

' / 
y / yy 
',y y 

yy 



9 12 
Md(0) [M solar] 



15 



18 



FIG. 5: The energy vs. the initial disk mass(MD(0)) for various initial angular momentum 
parameters (a (0)). 
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5(0) 

FIG. 6: The upper limits of the energy out of the system as a funtion of a(0). 
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FIG. 7: The effect of the magnetic field on the total energy out of the system for Md{0) 
3M0 (solid line) and Md(0) = TM© (dotted line) and M£,(0) = lOM© (solid - dotted line). 
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